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^ . Abstract 

^ ; 

■ We evaluate the one-loop corrections to the Wilson coefficient of the kinetic operator 

Q , in the operator product expansion of the differential B — > Xciv' decay rate. With a 

0^ I moderate cut on the lepton energy, the one-loop terms change the kinetic operator con- 

' tributions to spectral moments by about 20%. This amounts to a small correction for 

. leptonic and hadronic moments, except for those which vanish at the lowest order, where 

^P. \ the effect can be sizable. Together with a two-loop calculation of the leading-power rate 

' and an evaluation of the one-loop corrections to the Wilson coefficient of the chromo- 

^ ■ magnetic operator, our results will allow for a high-precision determination of \Vcb\ and 

^ i the b- and c-quark masses. 



1 Introduction 



Inclusive B — s> Xciiy decays are a precise probe of the underlying b- to c-quark transition be- 
cause hadronisation effects are small and have a simple structure. These effects are suppressed 
by powers of the heavy-quark mass and given in terms of a small number of non-perturbative 
parameters. In the heavy-quark limit, the hadronic decay rate becomes equal to the partonic 
decay rate. The leading corrections are of order and are given in terms of two non- 

perturbative heavy-quark parameters, /i^ and /i^ which are the i?-meson matrix elements of 
the kinetic and chromo-magnetic operator respectively. Schematically, the decay rate takes 
the form [1-4] 

where p = ml/ ml. The coefficients /, g and k can be calculated in perturbation theory: 

/ = f^'\p) + - f'\p) + (-)' f^'\p) + ,etc. (2) 

The general structure of the expansion is the same for other observables, such as partial rates or 
moments of the decay spectrum, but the calculable coefficients /, g and k are different. For the 
total rate the kinetic corrections have the same coefficient as the leading order, k{p) = —f{p)- 
Also for other observables, such as partial rates and moments, the kinetic corrections can be 
obtained from the leading-power differential rate, but the relations are more complicated and 
are evaluated to 0{as) for the first time in this paper. 

To turn ([T]) into a precision determination of \Vcb\ one needs the values of m^, mc and the 
heavy-quark parameters. Since the same parameters enter moments of the decay spectrum, one 
can determine these parameters by measuring not only the rate, but also a number of moments. 
To this end, lepton energy moments and hadronic invariant mass and energy moments are 
measured [5-12]. Using the results of these measurements, several groups have performed fits 
of the theoretical expressions to the experimental data [13-16]. The theoretical expressions 
that are used in the fit include one- loop corrections to the leading-power coefficients f{p) 
[17-26] as well as the Poal-part of the two-loop corrections [26-30], while the coefficients 
g{p) and k{p) of the power corrections are known only at the tree level. In addition to the 
second-order power corrections proportional to /i^ and Pq, the fits also include the third- 
order power corrections, which involve two additional hadronic parameters, p^ and p^^. [31] 
(the fourth order corrections are now available as well [32]). This technique yields the most 
precise determination of \Vcb\ together with very precise determinations of the heavy-quark 
masses. Already now, the estimated theoretical uncertainties are somewhat larger than the 
experimental ones [15]. In the future the experimental uncertainty will decrease further: the 
BaBar moment measurements which were used in [14,15] were published in 2004 and are based 
on 50 fb"^ of data [5, 60, and the recently published Belle measurements on 140 fb ^ [11,12], 
but combined the two experiments have already collected more than 1 ab~^ of data. Also, 
based on the convergence of the perturbative series of the rate for r-decay and based on 

^Very recently, Babar has presented preliminary results for hadronic moments based on 210 fb~^ [33]. 
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the size of the two-loop contributions that arise when converting the theoretical expressions 
between different schemes used in the literature, it has been suggested that the theoretical 
uncertainties in the results of the moment fits might be underestimated [34]. Whether the 
uncertainties are reliable is an important question because the value of extracted from the 
fit is a crucial ingredient for the determination of \ Vub\ from inclusive decays. After imposing 
the severe cuts necessary to eliminate the charm background, the prediction for the B Xuiu 
rate behaves as with n ^ 10 — 15 [35,36]. The value and uncertainty of the extracted \ Vub\ 
are thus strongly correlated with the value and uncertainty of rrih. 

It is clearly desirable to increase the precision of the theoretical predictions. To achieve 
this goal, two ingredients are needed: the leading-power moments have to be evaluated to 
two-loop accuracy, and the coefficients of the power corrections proportional to fi^ and fiQ 
need to be evaluated to one loop. In this paper we take the first and simplest step in this 
direction by evaluating the coefficient of the kinetic operator to one-loop accuracy. Let us 
stress that, while it is demanding, also the 0{al) calculation of the leading-power moments is 
doable. A few years ago such a calculation looked prohibitively difficult, but in the meantime 
the necessary methods to perform it numerically have been developed [37, 38] . Indeed, the 
two-loop correction for muon decay fi Xf.vv, the QED equivalent of i? ^ X^lv, has 
been evaluated recently using this method [39]. We use the same numerical approach for 
our one-loop calculation, because the size of the expressions involved is such that an analytic 
calculation does not look feasible. Since the method is gauge invariant (we perform the 
calculation without introducing a gluon mass), it is also suited for the calculation of the 
corrections to the coefficient of the chromomagnetic operator. 

The kinetic corrections are obtained by expanding the leading-power 0{as) expressions up 
to second order in the small residual momentum of the 6-quark inside the i?-meson. To have 
a check of our results, we perform the calculation in two different ways. A straightforward 
and tedious way of obtaining the kinetic corrections is to expand the leading-power Feynman 
diagrams in the residual momentum before performing the loop and phase-space integrations. 
The resulting expressions are long and involve terms which are individually strongly infrared 
divergent. Another complication arises because the expansion produces not only standard 
phase-space integrals, but also derivatives of such phase-space integrals, which arise from 
cutting higher powers of propagators. A much more elegant and efficient way to perform 
the calculation is to expand the result for the leading-order differential rate in the residual 
momentum. In this way one obtains results for the moments in terms of integrals over the 
leading-power rate and its derivative. In fact, without experimental cuts one obtains simple 
algebraic relations between the kinetic corrections and the leading-power moments. As a 
byproduct of our analysis, we obtain the 0{as) corrections to the leading-power moments and 
reproduce the numerical results of [26]. 

As the two-loop corrections to the leading-power rate and the one-loop corrections to the 
/ig terms are not yet available, it is too early to perform a detailed phenomenological analysis. 
Instead, we present numerical values for a few reference values of the electron energy cut. For 
the moments which do not vanish, the corrections we calculate turn out to be small, below 
1% as long as the cut on the lepton energy is not too strong. For the moments of the partonic 
invariant mass (p^ — "ti^)", on the other hand, which vanish at the tree- level and leading power, 
the corrections are larger, of order 30%. We expect the corrections proportional to /i^ to be 
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more important than the kinetic corrections. For the tree-level rate, they are roughly a factor 
ten larger than the kinetic corrections. 

This paper is organized as follows. In the next section, we explain how the kinetic cor- 
rections are calculated using the operator product expansion (OPE). In Section [3] we give 
parameterizations for the phase-space and loop integrals which are needed to perform the 
calculation. These parameterizations map the integration onto the unit hypercube and are 
such that infrared divergences appear only in a single variable and can be isolated before 
performing the numerical integration. We present our numerical results in Section m In 
the same section, we give formulae for the conversion of partonic to hadronic moments. We 
also sketch discuss how to convert our results into different schemes, such as the kinetic [40], 
potential-subtracted [41], IS* [42] or the shape-function scheme [43]. 



2 Evaluation of the kinetic corrections 

In this section we briefly explain how we evaluate the B — ^ X^iiy decay rate using the OPE. 
The application of the OPE to inclusive B decays has been worked out quite some time ago 
and the reader interested in more details should consult the original references [1-4] or the 
textbook [45]. Our goal in this section is to outline the necessary steps to obtain the result 
and to discuss some of the technicalities which are encountered in the course of the one-loop 
calculation. The reader solely interested in the numerical results can skip to section m 
The B Xclv decay is mediated by the effective Hamiltonian 

Kff = ^Kfc J' Jl = ^VcbC^' (1 - 75) b Ii, (1 -^,)v. (3) 
The decay rate factors into a leptonic tensor L^,^ and a hadronic tensor W^^j, 

dF = ^^fW ^^^^^^ ^^^^^^ L^,{p,,p,)W'^^{pj,, q) , (4) 
where q = pi + Pu and d/i(p) denotes the phase space 

^^(^^ = i2nr^ 2E 

in = 4 — 2e dimensions. Since the differential rate is a finite quantity, we could set d = 4. 
However, individual contributions to the hadronic tensor contain ultra-violet (UV) as well as 
infrared (IR) divergences which we regulate by keeping e 7^ throughout. The spin-averaged 
leptonic tensor is 

L^^ = Tr [;^£7m(1 - 75) ^7i^(l - Ts)] • (6) 

The hadronic tensor is obtained by taking the imaginary part of the time-ordered products of 
currents 

W^^, = -2ImT^„ (7) 

where 

T,. = -ij d'qe-^^^^{Bips)\ T [jt(x) J.(0)] \B{ps)) . (8) 
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We work in the kinematics = Mbv^ and our states are canonically normalized. In analytic 
calculations the hadronic tensor is usually decomposed into five form factors, but we prefer 
to directly evaluate the relevant product W^yL^^" . Since we use dimensional regularization 
we need to specify how we treat 75 in d dimensions. A definition of the axial current in d 
dimensions suitable for our purposes has been given by Larin [44] and we adopt it for our 
calculation. 

The operator product appearing in ([H]) is expanded in a series of local operators which 
corresponds to an expansion of the rate in inverse powers of the 6-quark mass. To perform 
the expansion, one first removes a rapidly oscillating factor from the 6-quark field by writing 
it as h{x) = e~^^'"^hv{x) . Using heavy-quark effective theory (HQET) [45,46], all the matrix 
elements necessary to second order in the expansion can be reduced to 

(Okin) = ^^{B{pB)\KmX \B{pb)) = -^il, (9) 

{Or..,) = ^^{B{pB)\hJ-a,.G^''K \B{pb)) = fih ■ 

The HQET parameter /i^ is often denoted by — Ai and is not renormalized. Up to terms 
suppressed by three powers of the heavy-quark mass, the decay rate thus takes the form 
q2 \y 12 

dr = ^ ^ '''' dn{pe) d/i(p^) [C3{v,pi,py){0s) + Ckin(^^,P<?,Pi.)(C'kin) + C^s.giv,pi,py){0^^g)] . 

The Wilson coefficients Ci{v,pi,Pu) of the three operators are independent of the external 
states and can be calculated using partonic initial and final states. To extract the coefficient 
C-kin{v,pe,Pu) of the kinetic operator, it is simplest to use an on-shell 6-quark with momentum 
Pb = rribVn + r^, which amounts to calculating the partonic decay rate b Xciu. To find the 
coefficient of the operator Okin with two derivatives, we expand the partonic rate to second 
order in the residual momentum r^. The result takes the form 

^ppartonic _ ^ + 4 ^ + A,, ^ r'^r" + 0(r^) . (10) 

TUb ^b 

At the loop level, the question arises whether to expand the diagrams before or after the loop 
integration. In general, either choice is valid as long as one evaluates the loop corrections to 
the operator product and to the matrix elements of the local operators Oi in the same way. In 
our case the situation is especially simple: since we perform the matching calculation on-shell, 
the one-loop corrections to the HQET matrix elements of the operators Oi vanish and the 
loop integration commutes with the expansion in the residual momentum. To have a check of 
our results, we evaluate the corrections in both ways. 

We can further simplify the calculation by averaging over the direction of the transverse 
momentum = — v ■ rv^. The component parallel to t'^ is fixed by the on-shell condition 
2mbV ■ r = — r^. Taking the average we have 

^ppartonic = A-A^ + V 4 T^(^?A«- - ^A^^V) + 0{r^) . (11) 

2m^ d — 1 
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Figure 1: Tree-level contribution to the hadronic tensor. 



To obtain the hadronic rate, we first bring the leading-power partonic matrix element into the 
form (Q by rewriting 

{b{p,)\K K IK'Pt)) = —{b{pi>)\kAbv \b{Pb)) 



{h{p,)\KfK\h{p,)) 



2ml 



{b{pb)\K K HPb)) 



(12) 



We then replace the partonic matrix elements by the corresponding hadronic matrix elements: 



dr =A—-{BipB)\Ki'b,\B{pB)) 
ZMb 



+ 

--A- 



d 



2Mf 



■{B{pb)\K {iDfb,\B{pB)) 



2ml 



i9, 



+ ... . 



(13) 



The ellipsis denotes terms which are suppressed by m^"^ or proportional to The operator 
Ojnag has a vanishing 6-quark matrix element and its coefficient is therefore not determined 
by our matching calculation. 

To illustrate the structure of the result, we now calculate the kinetic corrections at tree 
level. The tree-level contribution to the hadronic tensor is shown in Figure [TJ Taking its 
imaginary part and contracting with the leptonic tensor, the partonic tree-level rate is found 
to be 



^■pparto 



32 Gl iKbl' d/i(p,) d/i(p,) ph ■ p, {pb -q)-pi (27r) 6{{pb - qf - ml) . (14) 



We now expand up to second order in the residual momentum r^, average over the ±-direction 
and replace the partonic by the hadronic matrix elements to obtain the result 



dr = UtiGI \V,b? d/i(p,) d/i(p,) 



/o 5 [pI - ml) + /i 5' [pI - + /2 5" [pI - m^) 



(15) 
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where Pc = mf,v — q and 



/o = rribV ■ p^pc ■ Pi + 



2ml 3 



5mftV • piV ■ pi, - 2mbPi ■ p^ 



fi 



1^1 



V ■ PcV ■ Pu {2mbV + 5pc) -Pe-Pc-Pi (Sm^y + 2pc) ■ p^ 



f2 = 2nlmbV ■ p^p^ ■ Pi [{v -Pcf - pI] . 

The one-loop results for the rate have a similar structure, also in this case the result contains 
_ ^j^]-^ J7, = 0, 1, 2. To calculate moments of the decay spectrum, we introduce an 
integration over the partonic phase space 



1 



dpi 
(27r) 



d/i(p^) {2ttY S'^irUhV - q - Px) 



The result for the decay rate then takes the form 

d 



[dn] /o 



dmi 



[dn] /i 



d{m. 



2^2 



[dn] /2 



(16) 



(17) 



where we have used the notation 



[dn] 



[dn 



d/i(Px) J dfi{pi) J dfj.{p^) {27iy6'^{mbV - p^ - Pi - Pu) 



It turns out that for the total rate the derivative terms in (jTTj) do not contribute. However, for 
partial rates or spectral moments these terms do give non-vanishing contributions. To evaluate 
(|T7I) numerically, we need a suitable parameterization for the phase-space integral ( ITSi) . The 
necessary parameterizations, both for the tree-level phase space and the phase space with the 
emission of an additional gluon, needed for the 0{as) corrections to the rate, are given in the 
next section. 

An alternative, more elegant and efficient way of evaluating the kinetic corrections was 
described in [45] . Instead of expanding the diagrams which contribute to the hadronic tensor, 
one takes the result for the differential partonic rate and expands in the residual momentum. 
To derive this result, it is convenient to introduce the dimensionless variables x = 2Ei^/mb, 
y = 2Ei/'mb and = q^ /ml. When expanding in the residual momentum, one has 



X 



X H r ■ p^ 

nib 



y ^y-\ r-pi 

rrib 



-2 -2 

q ^q 



(19) 



Expanding to second order and averaging over the ^-direction, the hadronic differential rate 
is equal to [45] 



dr 



dx dy dg^ 



1 + 



2ml 



d 



d 



ox oy 



1 2 92 1 2 52 

— X 1 — y 

3 dx"^ 3 dy"^ 



^■pparto 



I '^ixy 2f) Y 
3 dxdy J \ dxdy dg^ 



(20) 
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Figure 2: One-loop contributions to the hadronic tensor. 

The /i^-term without derivatives comes from expanding the matrix elements, see (1121) . What 
makes this result particularly useful is that we can explicitly evaluate the derivatives using 
integration by parts when calculating moments. For the moments with a cut on the lepton 
energy 

[^'^2/'"(^')%o = / dx dy fy^., X- mQ{y - yo) , (21) 

one finds 



2 y-l partonic 



+ A_ [((^ + rnf + 2m + 2n - 3) y"' {qj - Amnx^-^y^~^ ^.2y+i] partonic 



+ ^ [((^ + + l)xyo - W) x^-'y^ifY 5{y - y,) + x^'y^+\e)'5'{y - yo)]'''^'"" • 

(22) 

The terms in the third line are boundary terms and vanish when setting = 0. In this case 
the kinetic corrections to the moments follow via simple algebraic relations from the leading- 
power term. The explicit relations for the moments we are interested in are given in Appendix 
Rl In the general case with a cut on the lepton energy t/q 7^ 0, one needs to also evaluate 
moments of the partial rate and its first derivative. To evaluate these boundary terms, it is 
important to keep in mind that g^, x and y are not completely independent: the rate includes 

2 

a factor 6{xy — (f') and at tree level the variables fulfill 1 — ^ + (f = x + y. More generally, 

when choosing a phase-space parameterization to evaluate (!22|) . the variables if and x become 
functions of y and the derivative in the third line of (!22l) then acts not only on the rate but 
also on the factors x°'{cp'y . 



3 Phase-space and loop integrals 

We now derive phase-space representations which are well suited for the numerical calculation 
of the kinetic corrections at one loop. The diagrams contributing to the hadronic tensor are 
shown in Figure O Their imaginary part receives contributions from virtual corrections as 
well as real-emission contributions. In the real-emission contributions, the imaginary part of 
the diagram is generated by an intermediate state with an on-shell gluon and on-shell charm 
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quark, so we need a parameterization of the b — > c+g+i+u phase space. The loop integrations 
necessary to evaluate the virtual corrections contain UV as well as IR divergences. The real 
contributions are ultraviolet finite, but contain IR divergences which cancel against the IR 
divergences of the virtual corrections. Because the quarks are massive, soft gluons are the 
only source of infrared divergences at one loop. 

To allow for a simple numerical evaluation, we map the phase-space and loop integrations 
to the unit hypercube. Also, since we want to calculate the rate and spectral moments with 
a cut on the lepton energy, we choose a parameterization in which the lepton energy is one 
of the variables. A last requirement is that we want the infrared divergences to be restricted 
to a single variable, so that they are easily isolated. It is convenient to split the phase-space 
integral into a hadronic and leptonic part 

J [dn,_e+g+^+p] = 1^1 [dl^b^Mp] J [dU,^c+g] ■ (23) 

3.1 Three-body phase space b ^ x + £ + u 

We denote the phase-space integration variables by Aj G [0, 1] with i = 1 . . . 4. Neglecting the 
lepton masses, we choose the momenta as 

Pb= K, 0,0,0), pe = {E,,0,0,Ee), = {E,, E, sin 6^,0, E, cos 6^) , (24) 

and parameterize 

T? y Z7 (1 -p-y) (1 - ^2) /3 0\ 1 fO^\ 

Ei = mb-, E^ = mb — , cosa*! = 2A3 - 1 , (25) 



with 



2 ' 
mf 



2 

K = 1 - (1 -cosei)y/2. (26) 



In terms of these quantities, the phase-space integral in d = 4 — 2e reads 



2 



271 



2d+i(2vr)2rf-2 
with the (i-dimensional solid angle 



dy / dX,dXs{l-p-yy-'^K'^-'{y{l-X2)y-''{{l-X3)X3)-\ 
^0 

(27) 



In the presence of a cut on the lepton energy y > yo, the transformation y = {1 — p — yo)Xi + yo 
maps the integration to the unit cube. It is simple to obtain the tree-level phase space from 
(1271) . To this end, one multiplies with 2'7i6{pl — ml) = 27r6 ((1 — p — y)X2) and integrates over 
A2. 
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3.2 Two-body phase space x ^ c + g 

We split the gluon three-momentum into a part in the direction of Px and an orthogonal part 
Pg = {Eg, 0, 0, 0) + Eg COS 9^ (0, + Eg sin 02 (0, p^) (29) 

with p±-px = and = 1. Expressed in terms of these quantities, the two-body phase space 
is 

For a given angle 6'2, the gluon energy is 



^ 2(E,-cos^2|p;|) ■ ^ ^ 

Note that the denominators of the real-emission diagrams 

{Vh - Vgf -^1 = -'^rribEg , (32) 

{Pc + Pg f -ml=pl-m\ (33) 

are independent of p^. The only dependence on p^ arises from the scalar products with 
lepton momenta in the numerator of the diagrams. The integration over the unit vector p^ is 
therefore trivial. The only non-vanishing integrals we need are 

C?'-V|l,l>ll>i| = |l,^5*^jf^.-2, (34) 

where 5*-' is the metric on the (d — 2)-dimensional sub-space, with b\ = d—2. For the evaluation 
of the diagrams it is simplest to parameterize the vector pj_ as 

p^ = sin ^3 (0, 0,1,0) + cos ^3^(0, E^ cos 9i + Ei, 0, -E^ sin ^i) . (35) 

\Px\ 

The integrand is then a second-order polynomial in cos ^3 and sin ^3 and the integral over p± 
takes the form 

j rf'^-v|l,cos2^3| =f^d-3 y"'rfcos^3sin'^-'^3|l,cos2^3|=f^d-2|l,^|. (36) 

To calculate the rate, we combine (127|) with ( 130|) and (136|) and rewrite cos ^2 = 2A4 — 1. 

The point A2 = corresponds to the kinematic configuration where soft singularities occur, 
since Eg oc p^ — ml = {1— p — y)\2 — >■ 0. Both propagator denominators (132|) are proportional 
to A2 and vanish at this point. The phase space (l30l) itself is proportional to Xl~'^^ so that the 
infrared divergences take the form 

1 1 



Xl+'' 2e 



5(A2) + 



A, 



Oie) . (37) 



The above relation is easily implemented into the code for the numerical evaluation of the 
diagrams. We evaluate both the divergent and the finite part numerically and check that the 
1/e divergences cancel in the final result within our numerical accuracy. 
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3.3 Loop integrals 

The virtual corrections involve loop integrals 

fj T J \ - [ fldu {1, fc^, k^K} 

y fc2(2pfe-A; + P)(2p,-A; + F+p2_^2) y-"^) 

with p1 = ml. We need the loop integrals for 7^ because we replace — > + r'^ and 
then expand in the residual momentum r^. The Feynman parameterization of the integral has 
the form 

Jo 

!^-l,v{up^,+ up'^), -v" {upl + up^:){upl + Mp^) + (^^'^l (39) 

with u = 1 — u and 

A = v'^ [mlu^ + m^-u^ + 2pf, ■ pcuu] + vu (1 + -uf ) (m^ — p^) . (40) 

Since we expand around the mass shell, the integral over the Feynman parameter v can always 
be done analytically after expanding, because the v integration completely factors after setting 
p"^ = m^. For the scalar integral, the v integration produces 1/e infrared divergences. The 
u integration on the other hand is always finite and done numerically, together with the 
integration over the tree-level phase space. 

4 Results for the moments of the differential rate 

When doing the calculation, it is simplest to evaluate moments using partonic variables. To 
distinguish partonic and hadronic quantities, we denote the partonic energy and invariant mass 
by and p"^, while writing Ex and p\ in the hadronic case. For the tree-level diagrams E^ is 
the energy of the charm quark, while in the diagrams where a gluon is emitted E^ = Ec + Eg. 
The partonic moments for which we present results in tables are defined as 

[w{EuE,,pI)\ = dEi jdE^dplj^j^^w{EuE,,pl). (41) 

We consider lepton energy moments w = E^ = {Ei/rrih)'^ with n = 1 . . .3 and the partonic 
energy and invariant mass moments w = -E" {p^ — p)™ with n + m < 3 and p = ml /ml. Note 
that we do not normalize the partonic moments to the rate. 

Numerical results for the moments without a cut and with a cut Ei > 1/4.6 (corresponding 
to El > 1.0 GeV for = 4.6 GeV) are shown in tables [T] and [2] for ^/p = mc/mt = 1/4. In 
the last column, we indicate the relative size of the kinetic 0{as) corrections. To estimate 
the relative size we use as = 0.22 and p^ = 0.4 GeV^. As explained in Section [21 we calculate 
the moments with two different methods and check that the results agree within numerical 
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1 


7T 




rv U?' 


% 


1 


0.6319(4) 


-1.123(4) 


-0.6319(6) 


1.125(8) 


0.1 


El 


0.1941(1) 


-0.348(1) 


0.0000(3) 


0.000(3) 


0. 


Ef 


0.06509(5) 


-0.1186(5) 


0.1085(1) 


-0.198(1) 


-0.2 


EJ 


0.02308(2) 


-0.0429(2) 


0.09232(5) 


-0.1714(7) 


-0.5 


Ex 


0.2667(2) 


-0.454(2) 


-0.6319(2) 


1.124(3) 


0.3 


El 


0.11576(9) 


-0.1845(9) 


-0.3667(1) 


0.610(2) 


0.4 


El 


0.05148(4) 


-0.0744(4) 


-0.17468(6) 


0.2534(8) 


0.4 


iPl - P) 





0.05693(3) 


-0.7305(2) 


1.281(3) 


-41. 


(pI - pf 





0.005754(3) 


0.20337(5) 


-0.5712(9) 


-19.4 


iPl - P? 





0.0011438(6) 





0.036918(7) 


23.4 


ExiPl-p) 





0.02970(2) 


-0.20013(6) 


0.2544(8) 


47.2 


EM-P? 





0.003373(2) 


0.09285(2) 


-0.2455(4) 


-17.1 


EM-p) 





0.015856(8) 


-0.03570(2) 


-0.0208(3) 


-1.8 



Table 1: Coefficients of the perturbative and power corrections to the moments (14T|) without a 
cut on the lepton energy for rric/mb = 1/4. Perturbative corrections are given in units of oisj-n^ 
the power corrections in units of yU.^/(2m^). All entries need to be multiplied by the common 
factor V"cfepm^/(1927r'^). The numbers in the table correspond to the partonic moments in 
the pole scheme. The last column gives the relative size of the kinetic 0[a^ corrections for 
default values of the parameters, see text. 

precision. In the tables, we also include numerical results for the tree-level moments. Using 
(1221) it would be simple to evaluate them analytically. 

The one-loop kinetic corrections are small for the moments which get contributions at 
leading power but sizable for the moments of — p)"". For example, the moment of j)\ — p, 
gets a correction of —40%. Compared to the tree-level contributions of the kinetic operator, the 
one-loop terms are typically suppressed by a few times a^/ vr. We thus expect that the extracted 
value of p^ will be shifted by about ±20% by their presence. Given the size of the kinetic 
corrections to the rate, we do not expect that the Oia^ corrections will affect the extracted 
value of I Kfe]. In the fit of Ref. [15] the value of p^ is varied by ±20% to estimate the theoretical 
uncertainty. The corrections we calculate are indeed of this size, except that varying the value 
of p^ correlates the change in all moments, while the perturbative corrections are different in 
each case. In the ffi of [14], the corrections are underestimated to be ^^qcD/'^fe ~ 0.0002: 
the contributions we find are roughly ten times larger. 

In Table [3] we give the result for the one-loop kinetic corrections as a function of the cut 
energy Eq and the ratio p = ml/ ml. To this end, we perform a quadratic fit around default 
values mc/mb = 1/4 and Eq = mb/4: ^ 1.15 GeV. The accuracy of the quadratic ffi is a 
few per cent except in cases where the corrections become very small. Tables with precise 
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0.09269(5) 


-0.1722(7) 


-0.6 


Ex 


0.2111(1) 
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El 


0.08917(7) 
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0.01801(1) 
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0.0015307(10) 
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-0.1897(3) 


-20.1 


EliPl-p) 





0.009147(6) 


-0.05271(2) 


0.0304(3) 


12.4 



Table 2: Coefficients of the perturbative and power corrections to the the moments (14T]) with 
A.GEi > 1 and mc/mi, = 1/4. Perturbative corrections are given in units of as/vr, the power 
corrections in units of /i^/(2m^). All entries need to be multiplied by the common factor 
G^FlKfeP'^fe/(1927r^). The numbers in the table correspond to the partonic moments in the 
pole scheme. The last column gives the relative size of the kinetic 0{as) corrections for default 
values of the parameters, see text. 

numerical results for arbitrary cut energies and charm-mass values can be obtained from the 
authors. 

Instead of the partonic moments, experimental papers present results for the normalized 
hadronic moments 

f-') = L iE^i^iE, y'(EM.). (42) 

To translate the results to hadronic kinematics we note that leptonic quantities are identical 
on the hadronic and partonic level. Using that the 5-meson momentum is = Mbv^, it 
follows that 

Ex = MB-v-q = MB-m^ + Ex, (43) 
P]c = {PB - qf = pI + 2Ex{Mb - rrib) + [Mb - rribY . (44) 
With these two equations, it is straightforward to translate our partonic results into hadronic 
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Table 3: Dependence of the coefficient of of the moments fj^Tl) on the lepton-energy cut 

b 

and the charm-quark mass. We define v = Amc/rrib — 1 and w = {AEq — mi,)/rnb and expand 
the moments to second order in these variables. The expansion coefficients in the table were 
determined by performing a quadratic fit to the exact results in the range 0.2 < mc/nib < 0.3 
and 0.5 < 4.6 En < 1.6. 



language. For example, the prediction for the lowest moments are obtained from the relations 



(Ee) = j^^rrib 



{Ex) = j-^{mb 



+ {Mb - mb) [1]) , 



(45) 



{pI -Ml) = — (ml [Pi - p] + 2mb{MB - rrib) [E.,\ + {{Mb - rubf + {ml - Ml)) [1] j . 

The moment with unit weight function is the rate: [1] = T{Ei > Eq), see ( I^Tj) . While the 
moments of p^. — ml vanish at tree level in the heavy-quark limit, the hadronic moments 
pI — M|, are nonzero. Using the above conversion formulae together with the results in Table 
[2] we obtain for example 



{P'x - Ml) 
{(pI - Ml)') 



0.860 + 1.59- 



-32.3 + 1.96- 



71 J 2ml 



0.939 + 7.00— 

71 



117.4-178.2— 



71 J 2m^j 



GeV^, 
GeV^ 



(46) 
(47) 
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Finally, let us note that we performed our calculation in the pole scheme, but it is simple to 
convert our result into different schemes. The pole scheme is calculationally most convenient, 
but plagued by large higher-order corrections. The problem arises because the definition 
of the parameters in this scheme relies on on-shell quark states, a concept not meaningful 
beyond perturbation theory. The resulting bad perturbative behavior can be improved by 
using parameter definitions with less infrared sensitivity such as MS quark masses. This is 
appropriate for the charm quark which we treat as light, however, the MS mass definition is not 
suitable for the 6-quark, because it is not consistent with HQET power counting. A number 
of alternative schemes, appropriate for heavy-quark processes, are available: they include the 
kinetic [40], the potential-subtracted [41], the IS" [42] and the shape- function scheme [43]. To 
one-loop order, the scheme changes from the pole into the new schemes have the form 

2 

mh = mhifif) + fif—ci + —^—C2, /i^ = /i^(/i/) 1 + — C3 C4, (48) 

where coefficients c, depend on the scheme. For example, in the kinetic scheme Ci = ^Cp, 
C2 = Cp, C3 = 0, C4 = Cp and /i/ = 1 GeV. The exact choice of the factorization scale /x/ 
is a matter of convention. When performing the moment fit, it is worthwhile to check to 
what extent the fit results are independent of this choice. A compendium of two-loop scheme 
conversion formulae can be found in [47]. Let us note that only the kinetic and the shape- 
function scheme provide improved definitions for the parameter /i^. To obtain the scheme 
change to 0{as) one simply replaces the pole-scheme parameters by the redefined ones in our 
one- loop results. Additional 0{as) corrections are generated when performing the scheme 
change in the tree-level results. We refrain here from explicitly changing the scheme, but the 
default values of our parameters are chosen such that they correspond to values which are 
typical for the improved schemes. 



5 Summary and conclusion 

We have evaluated the one-loop perturbative corrections to the coefficient of the kinetic oper- 
ator in the operator product expansion of the decay B X^lv. The corrections are typically 
(1 — 3) X ^ ^ 10 — 30% times the leading kinetic power correction. We thus expect that 
these corrections will change the extracted value of /x^ from the moment fit by about 20%. 
Whether this in turn has an effect on the extracted mf, and values is hard to estimate 
without performing the global fit. Since the kinetic corrections are very small for the total 
rate, the value of \Vch\ will likely not be affected. 

With the same numerical methods used here the one-loop chromo-magnetic and the two- 
loop leading-power corrections can be calculated as well. Once these are known, the theoretical 
precision on the predictions for the B Xciu decay will be superior to the experimental 
accuracy. These results will increase the precision of the extracted parameters and provide a 
nontrivial consistency check on the experimental data used in the fit. The calculation will also 
answer the question whether the currently used theoretical error estimates on the extracted 
parameters are realistic. This is of particular importance because the value and uncertainty 
of rub is a crucial input in the determination of \Vub\ from the inclusive B — ^ Xuiu decay. 
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A Moment relations 

Without cuts on the available phase space, the kinetic corrections to the moments are directly 
related to the leading-power moments. The relations can be derived from the general result 
(1221) . For convenience, we list here the explicit form of the relations for the moments we are 
interested in. We write the relations as A = B which should be read as 

r jppartonic 

For the energy moments the relations read 



El = + ^ (-^P' + - 2E,m^ , El ^ El + {AEI - 3Elm, - 2plE,) , 
and for the partonic invariant mass moments the relations are 
pI=pI + ^^ i-Pl + '^^-'^b - 2ml) , 

{Plf - {Plf + ^ Hpir + QEMfm, - 14(p^)^m,^ + ^Elplml) , 



E^pi = E^pI + ^ (^~plm, + yS>fe - 2E,mi ) , (51) 



EMr - EM? + ^ (-y (P^^-, + ^Elplm, - ^E^plml + 



EIpI ^ EIpI + ^ (-^(Px)' + \eIpI - ^E^plm, + ^Elm, - 2Elml 
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